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摘要 ”本文 给 出 一 元 半 正 定 多 项 式 一 种 降 次 平方 和 表示 方法 , 并 给 出 从 已 知 半 正 定 多 项 式 得 到 其 降 次 平方 
和 表示 的 算法 . 在 第 四 节 , 我 们 把 这 种 “ 降 次 平方 和 ”表示 思路 和 算法 用 到 多 元 多 项 式 上 也 取得 成 功 . 
关键 词 ” 半 正定 多 项 式 , 平方 和 表示 , SOS, 构造 性 , 希 尔 伯 特 第 17 问 题 


1 引言 

除非 特别 指明 , 本 文 讨论 的 多 项 式 和 算 阵 都 假定 是 实数 域 R 上 的 . 

设 p(z) 为 一 元 实 系数 多 项 式 . 如 果 p(x) > 0,Yz € R , Wp) 为 半 正 定 多 项 式 ， 如 
果 p(x) > 0, Vx € R, 则 称 p(x) 为 正定 多 项 式 . 将 半 正 定 的 一 元 多 项 式 表 示 成 一 些 多 项 式 
的 平方 和 和 (SOS, Sum Of Squares of Polynomials) , 是 Hilbert 第 17 问 题 的 一 部 分 . 容易 
证 明 (参见 [1]〉 ple) 半 正 定 的 充分 必要 条 件 是 存在 实 和 矩阵 五 使 得 p(z) = X(HHT)XT, 其 
中 X = (z"axt0i..9.,2,1) WB = 万 万 T 被 称 作 p(zZ) 的 Gram 矩阵. [1]rrol ETE SE Gram^4B 
阵 B 给 出 了 p(x) 的 9S0S 的 一 种 算法 . 
在 本 文中 我 们 证 明 互 可 以 取 成 下 三 角 和 矩阵 (甚至 在 很 多 情形 可 以 取 “ 空 心 的 下 三 角 矩 
WE") , 并 通过 直接 计算 及 (而 不 是 B) 给 出 一 种 SOS 算 法 . 在 该 算法 里 , 只 需要 解 一 个 自然 
三 角 化 了 的 线性 方程 组 , 所 涉及 的 运算 均 为 有 理 运算 , 只 涉及 很 初等 的 理论 知识 , 大 学 低 年 
级 的 学 生 甚至 中 学 生 都 能 理解 . 与 [1] 比较 , 本 文 算法 不 需要 计算 矩阵 的 特征 多 项 式 , 不 需 
要 做 实 对 称 矩 阵 对 角 化 , 也 不 需要 使 用 关于 多 项 式 根 的 Descarte 符 号 法 则 . 在 做 半 单 代数 集 
计算 时 , 我 们 的 方法 只 需 计 算 一 个 多 项 式 《有理 函数 ), [1] 中 需要 处 理 多 个 多 项 式 相关 的 
问题 . 

关于 一 元 半 正 定 多 项 式 的 其 它 SOS 表 示 , 可 以 参阅 [3]. 关于 多 元 (包括 一 元 FEE 
多 项 式 的 SOS 表示 的 较 新 进展 的 综述 性 论文 , 可 以 参阅 [5, 6, 10]. 
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B1005104) 资助 项 目 
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2 “理论 与 算法 

定理 1: (存在 性 定理 ) 一 元 多 项 式 p(x) =Y g cnx* 半 正定 的 充分 必要 条 件 是 它 能 表 
示 成 p. ， 

p(z) - 7 (Y 2) (1) 
i=0 \j=0 

其 中 Qj; >0,1<i<n. 
WEBH: 充分 性 显然 成 立 , 因此 只 证 明 必 要 性 . 

不 妨 假定 p(x) 首 项 系数 为 1, 由 半 正 定性 , p(z) AnKa E bpi, k — 1,2, mn, 
其 中 ax, Bk 均 为 实数 〈 半 正定 多 项 式 的 实 根 一 定 是 偶数 重 的 , 它 的 任意 一 个 2 重 实 根 我 们 可 
以 看 成 是 这 里 6 = 0 时 的 特殊 情况 ), i 是 虚数 单位 , 满足 这 = —1, 从 而 


TL 


p(z) = [T [e - ar)? +88] = sx (2) + ao). (2) 


k=0 


其 中 sn(z) = [kolt 一 Qj) 是 一 个 次 数 为 n 的 多 项 式 , 而 
=D (am) + E (eaea) + + (ija). (3) 


1<k,l<n,k#l 
当 p(z) 的 根 全 为 实数 时 , qe EEF GEKO , 否则 g(x) 是 次 数 为 2n — 2 的 半 正 定 多 项 
式 (至少 一 个 Bi 不 为 0) . 当 p(zx) 正 定时 (所 有 的 Bi +0) , g(x) 也 是 正定 的 
上 述 已 证 得 的 事实 , 利用 数学 归纳 法 , 可 得 定理 的 必要 性 
用 和 矩阵 的 形式 来 看 定理 1 会 更 直观 . 
推论 2: 存在 性 定理 , 矩阵 形式 ) 一 元 多 项 式 p(z) = Y, e Ck2* 半 正定 的 充分 必要 条 
件 是 它 能 表示 成 


q 


一 


plz) 三 和 (万 厅 T)XT (4) 
其 中 X = (x"”,X” 1,.…… ,XY,1), mHAXTIEÉAE,'€B4NMA 


Gm. n 
Qn,on—1  Qn—1,n—1 


H= M : Rg (5) 


Qn,0 Qn-1,0  ^'' Qo0,0 
其 中 aii 之 0,1 <i<n. 
根据 定理 1 的 证 明 过 程 , 我 们 可 以 得 到 万 的 以 下 性 质 . 
推论 3: 如 果 p(z) 是 正定 的 , 则 所 有 的 a;; > 0; 如 果 p(Z) 是 半 正 定 的 , 则 存在 一 个 加 ， 
ži > io 时 gj; > 0, 而 i < ioj € i 时 ay = 0. 
推论 3 表明 , 矩阵 有 的 主 对 角 线 上 的 元 素 全 部 都 是 非 负 实 数 . 如 果 这 些 对 角 线 上 元 素 不 
全 大 于 0, 一 定 存在 一 个 最 左 的 0, 在 这 一 列 及 其 右边 的 所 有 列 的 元 素 都 是 0, 而 在 这 一 列 的 


2 


左边 的 主 对 角 线 上 的 元 素 全 部 大 于 0. 这 说 明 平方 和 表示 中 , 可 以 “ 断 尾 巴 ”, 不 能 “ 掉 脑 
袋 ”. 同时 , 推论 3 还 表明 ， Ja 所 述 的 平方 和 的 次 数 是 严格 递 降 的 , 即 右 端 的 平方 项 的 次 
数 依次 是 相 邻 的 偶数 2m, 2n 一 2,… ,2no, 其 中 no 是 满足 0 < no < n 的 整数 . 因此 (1) 不 会 出 
现 2z2 = a? 十 z2 a , 也 不 会 出 现 如 zs + z2 这 样 的 表示 ARTD ， 
据 我 们 的 算法 后 者 的 SOS 可 以 表示 为 z5 +g? = (a? — 1x)" 十 24 十 37z2. 

为 讨论 方便 , 我 们 将 (5) 中 的 待定 参数 分 成 三 类 


L tel dn Gn, n—13* ^" 5 An, 0; 4n —1,05 7 ^ ^ 5 41,0; 0,0); 
D :=(aiil0 <i < n), (6) 


V :=(@ijll «4 <n, 0< j « 1). 


*E—3S LABEL 58 LAURISÉ m 十 1 行 的 那 2n 二 1 个 未 知 量 , 是 我 们 待 求解 的 参数 ; 第 二 类 万 是 位 
于 对 角 线 上 的 参数 (去 掉 两 端 ), 3 这 些 参 数 是 半 自 由 的 ; 第 三 类 V 是 余下 部 分 (位 于 三 角形 内 
部 ) , 这 些 参数 是 自由 的 . 在 下 文中 , D > 0 是 指 D 中 每 个 变 元 取 值 大 于 0 , 类 似 地 V = 0 是 
指 V 中 每 个 变 元 取 值 等 于 0 


Algorithm 1 正定 多 项 式 SOS 算 法 
输入 : 正定 多 项 式 p(Z) = Y rg cuz. 
输出 : 形 如 公式 (1) 的 SOS. 
1: 比较 (1) 两 端 Z 各 次 项 的 系数 , 得 出 方程 组 


2 
C2m = Gn 


C2n—1 = 2an nann_1 


Cn — 2an, nan,0 dee 


Cn—1 = 2ü4 1,n—10n—1,0 Sp ees 


€; —201101,0 7 
2 
Co 一 Q0.0 十 "P 


2: AE AR TED > 0 下 , OUT ACD) kA ZEB LT 89k 个 未 知 量 . 最 后 一 个 方程 不 必 
解 a0,0, 只 需要 解 出 4002 := a; o- 
3: 取 非 空 半 代数 集 5 := (A002(D,V) > 0 > 0) 中 一 个 点 Do, WW， 并 计算 Lo = 
L(Do, Vo). 将 Lo, Vo, Do 的 值 回 代 到 (1). 
4: return 形 如 公式 (1) 的 SOS. 


Algorithm 1} 28325 H 


特 
我 们 不 妨 将 

方程 组 (7) 
程 才 首 次 出 现 ; 


站 


已 


前 n 十 1 个 方程 求解 也 变 得 
需要 解 出 4002 := ai 0“《〈 而 不 是 aoo, E 
推论 4: 首 1 多 项 式 P(z) 


次 数 不 超 过 2 的 正定 多 项 式 , 表示 成 降 次 平方 和 时 , 因为 不 含有 待 
EF 方 和 表示 . 利用 连续 函 


的 降 次 了 
推论 5: 次 
半 正 定 而 非 1 


如 下 特点 : 
20 除开 首尾 两 个 方程 外 ， 
在 该 方程 的 “…. 
fta, >0,1<k<n). 
可 以 将 p(x) 的 首 项 系数 规范 化 为 1, 从 而 aw = 1, 这 样 可 以 规避 掉 这 个 开 方 运算 , 并 且 使 
更 简单 . 最 后 一 个 方程 , 看 上 去 必须 作 开 平方 运算 , 其 实 不 然 ， 


的 5 的 非 空 


性 


由 推 


论 3 得 到 保证 . 在 Algorithm 1 中 , 一 个 重要 的 


征 是 待 求解 的 变 元 集合 工 中 的 变 元 位 于 矩 
其 称 为 工 -算法 . 


阵 互 的 第 一 列 和 最 后 一 行 , ERER L” É, 


1) 已 经 


”部 分 不 会 出 现 


自然 三 角 化 , 即 工 中 的 第 k 个 未 知 量 在 (7) 的 第 k 个 方 
工 中 的 第 k 个 未 知 量 在 第 k 个 方程 中 是 线性 的 ， 
中 间 的 2n 一 1 个 是 线性 方程 容易 求解 (注意 条 


第 一 个 方程 


E 
ER 


关于 an 是 2 次 的 只 是 计算 算术 平方 根 , 但 通常 我 们 


mi 


H 


ZN 


代 时 也 直接 回 代 4002, 这 样 完 全 避 开 开 方 运算 ) . 


— 忆 科 ckzx 正 定 的 充 要 条 件 是 半 代数 集 


(A002(D, V) > 0, D» 0) z 0. 


数 的 保 号 性 
类 数 大 于 2 的 正定 多 项 式 的 降 次 平方 和 表示 有 无 穷 多 个 . 


定 的 参数 , 所 以 只 有 


E, 我 们 还 有 


EXER! 


=> 


BWiXSOSSUJ 


青 形 按 以 下 方法 处 理 . 


Algorithm 2 半 正 定 多 法 


: return 上 一 步 


: 作 无 平方 因子 分 解 p(z) = 
: 将 g(z) 校 Algorithm 1 写 出 SOS. 


输入 : 半 正 定 


多 项 式 D(Z) 


= E» eyz*. 


输出 : p(z) 的 SOS 


(f (2))* 


中 得 出 的 SOS. 


glz), 其 中 g(x) 无 平方 因子 因而 是 正定 的 . 


推论 6: 次 


数 大 于 2 的 半 正 定 多 项 式 只 有 唯一 的 降 次 平方 和 表示 的 充 要 条 件 是 它 的 根 全 


在 Algorithm 1 中 , A002 := ag oZ D, V RÆ IHA H K LAO02(D, V). 工 中 的 各 参数 


也 是 D,V 的 有 理 函 数 . 由 连续 
得 4002 > 0 的 待定 参数 为 有 理 数 , 这 档 


来 的 oi ;, j < i 都 是 


THA, 唯一 的 例 儿 


推论 7: 设 p(zZ) 是 首 1 的 有 理 系 数 正 


的 保 号 性 和 有 理 数 的 稠密 性 , 我 们 可 以 取 那 些 使 
fF 在 p(x) 是 首 1 多 项 式 的 情况 下 , 我 们 可 以 保证 求 出 
是 aoo 未 必 是 有 理 数 , 但 我 们 可 以 保证 4002 是 有 理 


二 


正定 


定 多 项 式 ， We 


=0 
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其 中 除 qaoo 外 《但 a30 是 有 理 数 ) ， 其余 所 有 系数 aij 均 为 有 理 数 . 
推论 8: 设 p(z) 是 有 理 系数 半 正 定 多 项 式 , 则 存在 有 理 系 数 的 加 权 降 次 平方 和 表示 


; 2 
Md (Y 2) 
1 一 0 j=0 
所 有 系数 di(> 0),a,; 均 为 有 理 数 . 
3 ”计算 实例 


fl: 对 于 x eR, 求证 z6 — rö — 2x4 + r3 +r? +1 > 0. 


解 : 这 是 一 个 经 许多 数学 爱好 者 用 来 讨论 正定 性 判定 的 一 个 例子 [13], 以 下 的 SOS 表 
示 便 是 正定 性 的 一 个 “明证 ”. 


例 : 本 例 来 自 文献 [3]. 写 出 zs 十 273 一 18x? 一 12x 十 117 的 SOS 表 示 . 
解 : 表示 (9) 来 自 文献 [3]. 


y*405*—1R8x"-—12-- 117 


" 217046V? — 29107 89156\° — 6597622612523 (9) 
= I | x4 | : 
u 21315 21315 29107 13224160752075 


用 本 文 算法 得 到 一 个 表示 如 (10) 


z^ 4-22? — 182? — 12z 4-117 — (z? 4 x — 10) 4 (m 4- 4) 4 1. (10) 


15): 本 例 来 自 文献 [16], 写 出 z6 —22?--Ax' —62? 62? — Ax +2 的 SOS 表 示 . 
解 : 本 文 算法 得 到 一 个 表示 如 (11) 


2Z6 一 22z5 十 42z4 一 6z3 十 6z2 一 47 十 2 


= (P-a +a) (o r+ 2 | (: y | > (11) 


Algorithm 1 给 出 的 正定 多 项 式 平方 和 的 缺点 是 有 时 把 一 个 “简单 ”的 多 项 式 会 写成 
“复杂 ”的 形式 . 
例 : 写 出 zs 十 1 的 SOS 表 示 . 
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解 : 根据 算法 可 以 得 到 一 个 SOS 表 示 


j X 1 17\” 1 3807 
e+1= [ur 2 F 二 22 ] . 12 
m EN e 27" ^ 64 4^ t 4096 (12) 


fj: ( 半 正 定 多 项 式 ) 写 出 z6 — 62? --14z* — 18:7? 1-17 z? — 12r +4 的 SOS 表 示 . 
解 : 事实 上 


xê — 62? + 14z^ — 182? 173? — 12z 4A — (z? +1) (s - 1? (z - 2). 


RE 


半 正 定 的 , 按照 Algorithm 2 我 们 得 到 SOS 表 示 如 下 


zê -625 十 14z 一 18z3 二 17z2 — 12x + 4 = (a? — 32? 1 2x)” F (x? — 3x 1 2. (13) 


15): REEL RUE ez) A Er 一 x 十 1 的 SOS 表 示 . 

fg: 这 里 给 出 两 种 表示 , 表示 (1 和 是 预先 取 定 W = 0 给 出 的 稀 玻 表示 (空心 三 角形 矩阵 
表示 , V 中 的 10 个 变 元 都 预先 取 值 为 0, 节省 了 很 多 计算 量 ), 而 (15) 则 是 未 预先 取 定 自由 参 
数 的 值 的 非 稀 玻 表示 . 


gU 一 Z 十 1 
-3 
2 一 和 十 1 
(7 TE ia 1) «(e Dos 2n (15) 


但 这 种 类 型 的 “稀疏 ”表示 并 非 总 是 存在 . 
例 : rê — 2x5 -5z? — Ax -- RA SOS IG 


Dy 3 di /3 3y d 
3 2 2 
— | | ls 1 
(s 2 2+3) (« T (3 i) 8 (16) 


fa Ze “ZSEE” ASOS T. 


解 : 事实 上 , 假定 它 有 空心 三 角 和 矩阵 


形式 的 表示 , 则 按 Algorithm 1 算出 来 的 


(81* - 38€ + 144) s* 


A002 = (4 s? + (t? +16) së + (28? -- 88) sê + 
)s8 2251?))/(—64s?t?). 


+ (481* — 228? + 324) s? + (1615 — 120# 4 


其 中 s = a22,t . HFA48t* — 228 t? +324 > 0,1616 —120t* +225 ? = t? (40? — 15)? > 


0, 所 以 4002 < "o > 0,t > 0 时 恒 成 立 . 
对 于 次 数 比 较 高 的 多 项 式 的 SOS 表 示 , 主要 在 计算 
些 情况 下 做 一 些 特殊 处 理 , 可 以 减少 计算 量 . 


半 代 数 集中 的 点 时 耗费 很 多 时 间 , DE 


例 : 写 出 p(z) = 2x16 — 4x1? — 2x14 + Ax? 十 272 — r’ + 7x* — 93? — 7z? + 9x +6 
的 SOS 表 示 . 
解 : ”通过 求 p(x) 导 数 的 实 根 , 我 们 可 以 求 出 它 的 最 小 值 等 于 1, 由 此 容易 求 得 到 分 


解 p(zZ) — 12 (x? — x 4- 1)? (2x? — x + 5), 


其 中 2712 — r -- 546—276 2RTIALT HJIERE & IU, 


ACRI FH BU TEL AEYAGK H2? — x --5 的 平方 和 表示 , 然后 可 以 得 出 p(z) 的 SOS 表 示 如 下 . 
2416 x15 — 2314 + 4r! + 251? — 35 +7r — 93? — Tr? 97 +6 
3 1 1 1 1 1 1V 
—29 xê | 5 | 4 | 3 | 2 
bu "d OU UC NR S z) 
1 3 1 5 . 1V 
2 T 5 4 3 | 2 | 
(ze mosse z7 z” "un iei) 
1 3 1 3 9 7 iX 
2 6 x5 4 | 3 2 | 
(5: 2^ WU r — q0» x) 
1 1 1 1 [i^ 
2 (525 — za La? — ca? Lad 1 
s 2 Ej g^ eam 3 em 
9 I9 
2 x? 4 | 
(一 -3 16^ x) 
1 2 
2 
e- l, x) 
313 
1 
pa (7 - : 
pi; 
对 于 一 个 正定 的 16 次 多 项 式 p(z), 按照 Algorithm 1 直接 计算 , 会 出 现 120 个 待定 参数 , 假定 


每 个 参数 取 2 个 不 同 的 值 构成 一 个 网 格 进行 搜索 , 那么 网 格 的 点 数 会 达到 1.3 x 108, 


搜索 完 这 样 一 个 网 格 , 是 一 个 非常 艰难 的 任务 . 正定 多 项 式 一 定 是 可 以 取 到 最 小 值 的 ， 
; 则 p(z) — b 将 是 一 个 半 正 定 多 项 式 ， 


果 我 们 有 办 法 求 出 


2p(z) 的 最 小 值 b (必定 大 于 0) 


要 


如 
并 


且 p(z) — 5 将 有 偶数 重 的 实 根 . 设 p(z) — b = f(x)(g(z))?, 其 中 f 是 一 个 无 重 因子 的 正定 多 
项 式 , 它 的 次 数 至 少 比 p(z) 低 2 次 , 此 时 可 以 先 做 出 /的 平方 和 表示 , 然后 再 做 出 p(x) 的 平 


方 和 表示 , 往往 可 以 大 大 降低 计算 难度 . 这 个 例子 里 ， 


16 次 正定 多 项 式 ,， 当 我 们 知道 它 


的 


最 小 值 后 , 分 离 出 来 两 对 2 重 实 根 , 这 样 次 数 降 低 4, 变 成 一 个 实质 上 的 12 次 多 项 式 问 题 , 对 
应 得 待定 参数 减少 到 66 个 , 几乎 是 原来 的 一 半 . 而 同样 的 每 个 变 元 取 2 值 的 网 格格 点 数 只 
有 7.4 x 10?? 个 . 

不 过 正定 多 项 式 计算 最 小 值 是 很 难 的 , 这 种 方法 并 不 是 对 所 有 正定 多 项 式 都 有 效 . 


4 多 元 多 项 式 的 SOS 
对 于 多 元 多 项 式 的 S9OS 问 题 , 我 们 也 可 以 转换 成 一 元 多 项 式 进行 处 理 . 


às 


Algorithm 3 多 元 半 正 定 多 项 式 SOS 算 法 
输入 : 多 元 半 正 定 多 项 式 p(x,y,… ,2). 
输出 : plx, y, = ,2) 的 SOS 表 示 , 或 “不 能 表示 成 多 项 式 的 平方 和 ”. 
1: 按 如 下 方法 将 p(x,Yy,… ,Z) 转 换 成 一 元 多 项 式 g(t): Wn, = 1, jT X p(t"s,y,--- ,Zz) 中 t 
的 次 数 记 为 n, 一 1,…，, EX p(t"*, tv... iz) 中 t 的 次 数 记 为 ns 一 1. 4 


q(t) p res T dd" (18) 


2: dE E p(m, y, ,zz) 平 方 和 中 出 现 的 单项 式 列 表 , 并 计算 将 在 q(t) 的 平方 和 表示 中 出 
现 的 t 的 各 次 老 ( 按 降 需 排列 JU = (tho these Lt), 
3: Bitat) EA Pk) SOSH A: 


2X (Zost) . (19) 


4: Ve q(£) 0 P8 AN GG CPUS SCR ul Em. RADA CRK) 解 ， 则 返 
回 “不 能 表示 成 多 项 式 的 平方 和 ”. 如 果 有 和解 ， 则 将 解 代 回 上 式 得 到 g(t) 的 SOS 表示 
如 (19). 

5: 在 (19) 的 两 端 顺 次 施行 如 下 运算 


rem( -t= — z, t), =- , rem( -,t"» — y, t), rem( - ,t"* — x,t). (20) 


反 过 来 得 出 p(z,y…… 2) SOS Xm. 
6: return 上 一 步 中 得 出 的 SOS. 


注 : 对 于 在 (19) 中 SOS 仍 然 取 下 三 角 和 矩阵 形式 , 其 存在 性 可 以 通过 半 正 定 和 矩阵 
的 Cholesky 分 解 来 实现 [15]. Cholesky 分 解 是 说 , 对 于 一 个 半 正 定 矩 阵 B 一 定 存 在 下 三 
角 和 矩 阵 且 使 得 B = HHT, 并 且 玖 的 对 角 线 上 的 元 素 均 非 负 . WRB 正定 , 则 五 的 对 角 线 上 
的 元 素 均 大 于 0. 实际 上 ， 定 理 1、 推 论 2 和 推论 3 也 可 以 利用 Cholesky 分 解 得 到 . 

例 : (来 自 [1]) 将 3 元 多 项 式 f (x,y,z) = x8 +4ry?z + y? -2y*z? - y?z^ -- A28 给 出 
平方 和 表示 . 
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解 : 注意 z,y > 的 最 高 次 数 分 别 为 6, 令 


pae 


g(t) eS I 一 41288 | 1136 | 2114 | 4 199 


1 jz,y 纪 的 平方 和 表示 中 只 能 出 现 如 下 的 项 : x3, y3, y? z, yz? 22]. 从 而 9 的 的 平方 和 
表示 中 只 能 出 现 以 下 项 [#3,t21,457,t93,t129]. 设 g(t) 的 SOS 有 待定 形式 


rr 
FE 
| Gk 
H 


2 
g(t) = (a4t? + ast?” + aot + a t? + agt?) 


比较 g(t) 的 两 种 表示 的 关于 t 的 各 次 暴 的 系数 , 得 到 一 个 关于 这 些 待定 系数 的 方程 组 , 求 得 
方程 组 的 一 个 解 as 2, a3 0, a2 0,a4 0, o 0, bs 1, b2 0, bı =l, bo = 0, C9 = 
2,c1 = 0, co = 1, dı = 0, do = 0, eo = 0. 因此 得 到 g(t) 一 个 SOS 表 示 


44258 | t186 | 21114 | 4190 .. 142 .. 16 — (21129)? 4. (193 p | sd | t (21) 


在 上 式 两 端 顺 次 施行 如 下 求 余 运 算 


rem( - t — z,t), rem( - t — y, t), rem( - ,t — x,t). (22) 
则 得 到 原 多 元 多 项 式 f (zx,y,z) 如 下 SOS 表 示 


28 十 473g2z +y? 2y'z? t 42 = (222? c zy) (a? -2y?zy. (23) 


下 一 个 例子 也 是 来 自 [1], 原文 中 对 本 问题 的 结果 是 否定 的 , 即 判断 不 能 表示 成 SOS. 有 
关 不 能 表示 成 多 项 式 的 平方 和 的 问题 可 以 参阅 [2]. 

f5: (来 自 [1]) AH f(z,y,z) 2 xt + a?z -2oz7y? + z^ 的 90S 表 示 . 

解 : 因 z 的 最 高 次 数 为 4  z,y 的 混合 最 高 次 数 项 为 z2y2，z 的 最 高 次 数 为 4. ^ g(t) := 
f (t, 00,03) = 0? 08-21? e. 在 f 的 SOS 表示 中 会 出 现 的 项 为 [z2,yz,zz, 22]. 在 g(t) 的 SOS 
表示 中 会 出 现 的 项 为 [如 ,#5,t14,t29]. 假设 g(t) 具 有 以 下 形式 


g(t) — (ast?* -- aat! 4- a1t* -- agt?)* 


BE epu d" 


( 
z5 (cit? 十 ap^) 
- (dot?)*. 
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比较 g(t) 两 个 表达 式 的 系数 得 到 一 个 方程 组 , 求 得 方程 组 的 一 个 解 os = 1,a2 = 0,ai 
0, ao -= —1/2, b» IN 0, bo 1/2, cı — V2, Co = 0, do E 2/2. 从 而 我 们 得 到 SOS 表 示 


1 S 1 lr 1 
t9? 4. 18 9 t2 4 14 = (e ze) | (e | ze) +2t +3t. (24) 


在 上 式 两 端 顺 次 施行 如 下 求 余 运算 


rem( - , t? — z,t), rem( ,t? — y, t), rem(:,t— x,t). (25) 


JU f 58 JE E 76 e DURS f (or, y, 2) ll FSOSSEZS 


iN iN 1 
2 十 23z 十 27272 十 太一 G a | (= | z7) - 2x?y? 4 Tu (26) 


Bi OCKH[AD 证 明 f(x,y) ^ 2? c 22*y - 5a?y* c Azy? +y’ 7 0. 
解 : 因 z,y 的 最 高 次 数 都 是 6, 令 
g(t) = ft) m4 HAE BD 219 A tË. 


因为 在 f(x,y) 的 SOS 表 示 中 只 能 出 现 单项 式 [z3,z2y, zy2, y?] ,所 以 在 g(t) 的 SOS 表 示 中 只 能 
出 现 单项 式 [ 妇 , 妇 , 妇 5, 121). 假设 


g(t) = (ast?! + ast? + at? + a 

(bot!9 -- bt? 4 bot)? 
十 (cit? E cot?) 
+ (do£)? 
比较 g(t) 两 个 不 同 的 表示 的 系数 得 方程 组 , 求 得 方程 组 得 一 组 解 为 a3 = uas = 2a 
0, ao = —1/2, bz = 1,5 = 1/2, bo = —1/4, 1 = 3/2, co = 3/4, dj = 1/8. 所 以 我 们 得 到 g(t) 
的 表示 


t 十 4136 51:99 41912 4 16 


[2x 1 Pow 3 gc d 
= (elpa ;^) (e+e ze) ($e i6) uu (27) 


在 上 式 两 端 顺 次 施行 如 下 求 余 运算 


rem( - t — y,t), rem(:,t— x,t). (28) 
则 得 到 原 多 元 多 项 式 f(x,y) 如 下 SOS 表 示 
ToTY E TY 4 Amy? +y’ 


EIC 1 1 2^. /3 3 1 
二 3 2. 2 3 2.) 4 3| 1 6 (9 
(vs TY ze’) (zv ;* V iu) + (Boy 2 g7 (29) 


Bl: (Motzkin 的 著名 例子 ) 证 明 f(z,9y) = sty? + ryt- 3a2y? 十 1 不 能 表示 成 多 项 式 
的 平方 和 . 
解 : 由 算术 几何 平均 不 等 式 知 f(z,y) > 0 始终 成 立 . 


g := f(t, t) =t” +t a1 


h := [1, £, y, yx, a y ey] 


L:= [1,4 t,t, t,t] 


6 j 2 
gl:— y (x ——— 
j=1 \i=1 


比较 "和 91 的 系数 ， 我 们 容易 确定 出 a6 6 = 1, 6,5 = 0, 6,4 二 0, 6,3 = 0, 5.5 = 1, 5.4 = 


0, as,2 = 0, Q6,2 = -1/2 04,4? — 05, 一 3/2, Q6,1 三 —04,304,4, 05,3 = Ü, a43 = 0, a3,3 = 0. 此 
时 有 
91 一 9 三 (2a42a44 十 2 ds.1)t* + 2044044 P 
1 3 9 
T (1 asa“ RE 5 asa | 1 | 42.3" | 03,3" | ag E (30) 


+ (2 a2,102,2 + 2a3,103,2 + 2a4,104,2)t 


2 2 2 2 2 
+ 01,415 + 024^ o 034^ +04,1" + 054 1. 


E R30 A mte E] RAA BE S T0, 从 而 gl - 9 = 0 不 可 能 成 立 , Xi Harty? + a?y* 一 
32y? 十 1 不 可 能 表示 成 多 项 式 平方 和 . 


5 ”补充 说 明和 讨论 
oe 在 [7] 中 有 专门 的 计算 半 代 数 集 的 软件 包 SDPTools 可 以 供 大 家 使 用 , 另 可 参考 [7, 9). 至 
于 无 平方 因子 的 分 解 算法 , 则 可 以 参见 [8]. 
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